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Abstract 

We review recent theoretical developments about the role of spins, electron-electron 
interactions, and spin-orbit coupling in metal nanoparticles and semiconductor 
quantum dots. For a closed system, in the absence of spin-orbit coupling or of 
an external magnetic field, electron-electron interactions make it possible to have 
ground states with spin S > 1/2. We review here a theoretical analysis which makes 
predictions for the probability of finding various values of spin S for an irregular 
particle in the limit where the number of electron is large but finite. We also present 
results for the probability distribution of the spacing between successive groundstate 
energies in such a particle. 

In a metallic particle with strong spin-orbit interactions, for odd electron num- 
ber, the groundstate has a Kramers' degeneracy, which is split linearly by a weak 
applied magnetic field. The splitting may be characterized by an effective g-tensor 
whose principal axes and eigenvalues vary from one level to another. Recent calcula- 
tions have addressed the joint probability distribution, including the anisotropy, of 
the eigenvalues. The peculiar form of the spin-orbit coupling for a two-dimensional 
electron system in a GaAs heterostructure or quantum well leads to a strong sup- 
pression of spin-orbit effects when the electrons are confined in a small quantum 
dot. Spin-effects can be enhanced, however, in the presence of an applied magnetic 
field parallel to the layer, which may explain recent observations on fluctuations in 
the conductances through such dots. 

We also discuss possible explanations for the experimental observations by Davi- 
dovic and Tinkham of a multiplet splitting of the lowest resonance in the tunneling 
conductance through a gold nano-particle. 

*This article is slated to appear as a chapter in " Nano-Physics and Bio-Electronics", 
edited by T. Chakraborty, F. Peeters, and U. Sivan (to be published by Elsevier Co.). 
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1 Introduction 



During the last two decades the fabrication technology of small conducting 
islands, known as quantum dots, using semiconductor heterostructures, spat- 
tering of small metal grains and other methods have advanced so much that 
they can behave, under the right conditions, as artificial atoms [1] . In contrast 
to natural atoms, these artificial atoms do not have special symmetries, unless 
special experimental efforts are being performed[l]. 

The electron spin is responsible for a number of interesting effects in small 
chaotic conducting islands at low temperatures which are quite distinct from 
the role of spin in a bulk material. The role of spin is modified by electron- 
electron interactions in a way that has consequences for the distribution of the 
energy separations between ground states with different number of electrons, 
as well as for the probability of finding different spin quantum numbers at a 
fixed number of electrons. Other interesting effects are produced by spin-orbit 
coupling, which is important distribution of energy levels and wavefunctions 
in closed quantum dots and may modifies its properties, and which affects the 
conductance distribution for an ensemble of open quantum dots. The addition 
of source and drain leads, and in case of semiconductor heterostructures gates 
that control the charge on the dot, allows one to measure the properties of the 
dot-leads compound as a function of V, the difference between the potentials 
on the source and drain, and V g , the potential on the gate lead. 

In Sec. 2 of this paper we discuss a few effects of electron-electron interactions 
in a closed dot, with no spin-orbit coupling and no significant Zeeman field. 
The analysis reviewed in this section has been developed by a number of 
research groups in the last few years. We re-derive here an effective low energy 
Hamiltonian that was first discussed in Ref. [2] using renormalization group 
(RG) arguments. Then, we extend our studies [3] for the ground state spin 
of a quantum dot, and analyze its influence on the Coulomb blockade peak 
spacing that appear in the conductance as we sweep V g at zero bias voltage, 
V. Further details, concerning the parameters of the effective model in actual 
systems and the relation of the effective model to other models, are given in 
the Appendices to the paper. 

In Section 3, we give a brief review of recent work on the effects of spin-orbit 
coupling in metal nanoparticles and GaAs quantum dots. The reader will be 
referred to the literature for a fuller account. 

In Section 4, we discuss a problem motivated by experimental observations of 
Davidovic and Tinkham [4] of a multiplet splitting of the lowest resonance in 
the tunneling differential conductance through a gold nanoparticle as a func- 
tion of the source-drain potential, V. We consider two possible mechanisms 
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that may lead to effects of that type. The first is due to an almost degenerate 
ground state (Subsection 4.1) and the second due to a nonequilibrium phe- 
nomenon (Subsection 4.2). In Subsection 4.3 we compare these mechanisms, 
which both involve electron spin and/or electron-electron interaction in an 
essential way, and suggest experimental ways to distinguish between them. 



2 Electron-Electron Interactions in a Closed Dot 

When the shape of the dot is symmetric [1], a part of its single particle levels 
are degenerate, just like in an atom with a spherically symmetric potential. 
The exact electron many body configuration, which is set by the repulsive 
interaction between the electrons combined with the Pauli exclusion principle, 
is summarized in a set of rules, known as Hund's rules. The first of them 
assert that a partially filled set of degenerate levels will have the maximal 
spin that is consistent with the exclusion principle. This happens because 
for a system of several electrons the "most antisymmetric" coordinate wave 
function has the largest spin. The most antisymmetric wave function minimizes 
the Coulomb repulsion between the electrons because it vanishes when there 
are two electrons at the same point. Strong spin-orbit interactions may change 
Hund's rules. 

To create symmetric dots, special (experimental) effort is necessary. Any 
generic dot, however, does not have any special symmetry and can be con- 
sidered as a chaotic one, its single particle levels may assumed to be random, 
and are described by random matrix theory (RMT). Many theories have con- 
centrated on the statistical properties of random levels and the way they 
can be used to understand chaotic dots in actual experiments, for review see 
Ref . [5] . The interactions between electrons in the dot and between them and 
the environment were commonly described by the so-called constant interac- 
tion model. In this model all the effects of interaction are cast into a single 
capacitance that describes the change in the energy of the system due to the 
dot's charging. 

This class of models was very successful in explaining and predicting experi- 
mental results as Coulomb peak height fluctuations, conductance fluctuations 
through an open dot and so on [5]. However, it fails to explain experiments 
measuring motion of peaks in magnetic field [6], distributions of Coulomb 
blockade peak spacing [7-10], and multiplets that appear in a single Coulomb 
blockade peak [4] . Motivated by the failure of the constant interaction model, 
a few models that extend it were suggested [2,3,11-13]. In Subsection 2.1 we 
derive, by integrating out fast degrees of freedom in the RG sense, an effec- 
tive low-energy-Hamiltonian, Ti e s- This effective Hamiltonian, first discussed 
in Ref. [2], describes the properties of the quantum dot at energies smaller 
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than the Thouless energy of the system, which is inversely proportional 
to the time it takes to cross the dot along its largest dimension. We then 
show that a proper choice of parameters for 7i eff and its analysis reproduced 
other models for the effects of electron-electron interaction in quantum dots. 
Appendix A presents some details of the relation and equivalence of different 
models. Throughout this article we will neglect fluctuations of the interaction 
parameters, we discuss the limits of this approximation in Subsection 2.1.1. 

In the process of averaging over the fast motion of the electrons in the Fermi 
sea, three channels of interaction appear: The direct /charging channel, the 
exchange/spin channel and the Cooper channel. The actual values of the pa- 
rameters that describe these channels in 7i e ^ depend on the fast motion of 
electrons. We can have a richer situation when there are intermediate scales 
between the Fermi energy and the Thouless energy, e.g.: the inverse of the 
mean free time, r, between elastic collisions with static impurities in diffusive 
systems. In Appendix B we find, using the Fermi-liquid-theory, the effective 
exchange interaction parameter for several three and two dimensional sys- 
tems. Its dependance on electron density in two dimensional Si-MOSFET and 
AlGaAs ballistic heterostructures is estimated. Appendix C deals with the 
renormalization of the Cooper channel. 

After a detailed discussion of the model in Subsection 2.1, we use it to analyze 
several physical quantities. In Subsection 2.2 we present an extended version 
of our study in Ref. [3] of the ground state spin configuration distribution. The 
following Subsection (Subsection 2.3) discusses the effects of spin fluctuations 
in the ground state on the Coulomb blockade peak distribution, (see also 
Ref. [12]). In Subsection 2.3.1 we compare some predictions of the theory with 
available experimental results [9,10]. We find that, although several features 
of the theory are observed in experiments, there is still disagreement between 
theory and experiment. We will discuss the possible causes for this discrepancy. 

2.1 Effective Hamiltonian 

To find properties, such as the ground state spin distribution of electrons in the 
dot in the presence of electron-electron interactions, we would like to describe 
them at low energy in a simple form. Ref. [2] shows that in the limit of large 
dots the effective Hamiltonian, at energies smaller than -Erh, i s: 

H eS = ]T e^l^s + J S S-S + J C T^T + U C N 2 . (1) 

fJ,,S 

Here, N = E^ s =u V{ s Vw> S = fV^As'Vv and T = E^t^i- 

The ji sum runs over M = g = 27r£' Th /A states. The symbol A = \/{uV) 
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denotes the average single particle level spacing in a dot of volume V and 
thermodynamic density of states v (per spin). Notice that in our notation A 
refers to the level spacing for a single spin state in the dot. In a dirty dot of 
length L and diffusive constant D, E Th = D/L 2 , while for a single dot we 
replace D by ~ vfL. 

The first term in (1) is universal: is a single electron eigenenergy of a random 
matrix, and the operator ^ s is the creation operator of an electron with spin S 
at an eigenstate \x of a random matrix. The set {e^} depends on the symmetry 
of the problem. In the presence of time reversal symmetry the eigenenergies 
are taken from the Gaussian orthogonal ensemble (GOE, (3 = 1), and in the 
absence of time reversal symmetry from the Gaussian unitary ensemble (GUE, 
(3 = 2). When strong spin-orbit interactions are present they are taken from 
the Gaussian symplectic ensemble (GSE, (3 = 4)Q 

The direct Coulomb interaction constant U c , the exchange constant J s , and 
the interaction in the Cooper channel J c depend on the specific system and 
the model one uses for the interaction. In addition, there are non universal 
corrections to Hamiltonian (1) that vanish, however, in the limit of g — > oo. 
(See also in Subsection 2.1.1.) When the time reversal symmetry is broken the 
interaction in the Cooper channel vanishes. We will see below that even in the 
presence of time reversal symmetry the interaction in the Cooper channel is 
reduced due to a "screening" by fast electrons(see also Appendix C). 

To understand what are the effective low energy interaction constants U c , J s , 
and J c it is useful to describe the derivation of the effective Hamiltonian (1) 
in terms of a RG scheme. When the temperature decreases we integrate out 
progressively the fast motion of the electrons with energy far from the Fermi 
energy and find effective coupling constants in the direct (U c ), exchange (J s ) 
and Cooper (J c ) channels [14,15]. First, the fast motion of the electrons, at 
energies of O(E-p) away from the Fermi level, "dresses" the bare electrons and 
forms quasi-particles. The Landau Fermi-liquid theory describes this dressing 
process[14,15] and the way it renormalizes the system parameters. In Ap- 
pendix B we use this theory to estimate J s for several situations. 

The Fermi-liquid "dressing" continues up to energies of order 1/r. Below 1/r 
the motion of the electrons becomes diffusive and new diffusion singularities 
appear. In a situation where one of the dimensions of the system is much 
smaller than the others, the system may become quasi-two-dimensional at 
frequencies smaller than the Thouless energy that is related to the short di- 

1 The symmetry index (3 counts the degrees of freedom of the matrix elements 
of the single particle Hamiltonian, (3 = 1,2, or 4 if its elements are real complex 
or real quaternion numbers, respectively. A magnetic flux ~ 4>o/^ through the 
dot, where <po = hc/e ~ 4.12 x 10~ 13 TesZa • m 2 , leads to a time reversal symmetry 
breaking. 
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mension. This reduction in the dimension of the system enhances the diffusive 
singularities and changes the flow of the interaction parameters [16]. It appears 
that the RG flow in the Cooper channel is sensitive to disorder more than the 
RG flow in the other channels[16], (see also Appendix C). However, in cer- 
tain situations, especially when disorder is strong, we also have to consider its 
effect on the flow in the other channels[16]. 

Finally, we arrive at temperatures T < E^h and are left with an effective 
"zero dimensional" Hamiltonian, H e g- The length scale associated with such 
low temperatures is larger than the system size and therefore the interaction 
parameters are constants that do not depend on the site or state index. We 
note that in ballistic samples, i.e., when the electrons cross the sample before 
suffering substantial scattering from impurities we may use the Fermi-liquid 
theory, without additional complications due to the diffusive motion at inter- 
mediate scales. 



2.1.1 Fluctuations in the interaction parameters 

In practice, not all the samples have exactly the same shape and/or impu- 
rity configuration, we therefore expect to find sample to sample (mesoscopic) 
fluctuations in the RG process that will lead to fluctuations in the interaction 
constants of the different channel. By assumption (that is motivated and sup- 
ported by numerical analysis [17,18]) the sample to sample fluctuations in the 
single particle levels are described by random matrix theory. The eigenenergies 
and eigenfunctions are those of a random matrix. 

The random electron states have random charge distributions and hence we 
expect that the interaction parameters J c , U c and J s themselves will fluctuate 
with the electron number. The interaction in the Cooper channel is reduced 
by a logarithmic factor (see Appendix C) and we will neglect its fluctuations. 

When an electron is added to the system charge flows to the dot edges and 
leads to fluctuations in the self-consistent nonuniform potential [19]. These 
fluctuations scale as ~ (r s / ■ s /g)A, and give the largest contributions to the 
fluctuations in U c . [See Eqs. (B.3) and (B.7) for a precise definition of r s in two 
and three dimensions.] Fluctuations in U c are relevant to the Coulomb peak 
spacing distribution (Subsection 2.3.1) and to nonequilibrium effects (Subsec- 
tion 4.2) in the conductance through the dot at finite bias voltage. 

The short range part of the Coulomb interaction determines the exchange in- 
tegral in the expression for the exchange interaction parameter J s . We can 
therefore use a contact interaction model to find its fluctuations. The fluctua- 
tions in the interaction parameter in that case, for three dimensional samples, 



are [20] y/var J s xs J s max I l/(2y/2N), cs/g \ , where C3 is a numerical number 
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of order 1, and N ~ Ep/A is the total number of electrons in the dot. The 
first term is larger for g 2 ^> N ^ L < (kpl)l, with L the sample size, / the 
mean free path, and hp the Fermi momentum. For ballistic systems I ~ L, 
and the first term is much larger. Semiclassically [21], the first contribution 
arises from direct trajectories between two points in the dot, and the second 
is built from indirect trajectories with possible scattering on the surface or 
on static impurities. In two dimensional samples a similar calculation [12,13] 
gives, y/var J s pa J s max |c2 log N/ y/N, c' 2 /g} , where C2 and c' 2 are of order 1. 
For ballistic systems we find that, as in three dimensional samples, the first 
term is larger. 

The actual contribution of the fluctuations in J s to the fluctuations in the 
ground state energy is larger by a factor ~ y/g, as in the calculation of the 
ground state energy we have to include interaction of a single spin with ~ g 
electrons [12]. 

In any case, in the universal limit, when g, N — > oo we can neglect the fluctu- 
ations in the interaction parameters. 

2.1.2 A toy model with contact interaction 

In Ref. [3] we discussed a model with contact interaction, described by the 
Hamiltonian 

^toy = c n,s^o(n, m)c m>s + uM cl A ci d c nd c nA . (2) 

n,m,s n 

Now m runs over M coarse grain sites (in real space) and Ho(n, m) is a random 
matrix. The parameter u describes the strength of the interaction between the 
electrons. 

This toy model was analyzed within the self consistent Hartree-Fock approx- 
imation in the limit of large M [3]. We will show in Appendix A that to 
first order in u the toy model (2) is equivalent to the Hamiltonian (1) with 
parameters 

U c = u/4, J s = —u and J c = u. (3) 

Higher order corrections in u preserve the symmetry and the structure of 
Hamiltonian (1) but give different values for its parameters. In particular, 
as shown in Appendix C for positive u (which corresponds to a repulsive 
interaction) they reduce J c by a factor oc log M. 
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2.2 Ground state spin distribution 
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Fig. 1. The probability distribution P(S) of the ground state spin of a quantum dot, 
computed from Eq. (4) for different values of the interaction parameter A = — J s /A 
with time reversal symmetry. Solid histograms are for integer spins, dotted ones for 
half-integer spins. All the graphs are starting with spin increasing to the right in 
increments of 1/2 spin. 

In this section we will study the ground state spin of Hamiltonian (1) assum- 
ing that the interaction constant J c = [as it renormalized towards zero by 
electrons with energy larger than the Thouless energy and is further renormal- 
ized within the toy model (see Appendix C)]. In this case, the Hamiltonian 
(1) becomes noninteracting within each spin sector. For a fixed number of par- 
ticles A" we can neglect also the charging energy U c . The spin of the ground 
state is then found by minimizing the energy Eg{S) of the lowest lying state 
with total spin S, as a function of S. Since the lowest energy state with spin 
S has precisely 2S singly occupied states, all lower lying states being doubly 
occupied, one has [3] 

S-So 

Eq{S) — E G (S ) = ^ { £ N+n+2S ~ £n+i-h) 

+ J B [S(S+1)-S (S + 1)}. (4) 

where So = (1/2) if the total number of particles A^ is even (odd). Since 
the precise positions of the energy levels {e^} in Eq. (4) fluctuate from grain 
to grain, the ground state spin S does so as well. Plots of the probability 
distribution P(S) are shown in Fig. 1 and Fig. 2 for different values of the 
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Fig. 2. The probability distribution P(S) of the ground state spin of a quantum 
dot, computed from Eq. (4) for different values of the interaction parameter A 
without time reversal symmetry. Solid histograms are for integer spins, dotted ones 
for half-integer spins. All the graphs are starting with spin increasing to the right 
in increments of 1/2 spin. 



dimensionless parameter A = — J s /A. The effective parameter A, that includes 
renormalization from fast electrons, is calculated for small metallic dots and 
semiconductors dots in Appendix B. The distributions are obtained by taking 
the levels e M from the GOE {(3 = 1), or, when time-reversal symmetry is 
broken, from the GUE, ((3 = 2), and minimizing Eq. (4) with respect to S. The 
way spin-orbit coupling reduces the effect of the interaction in the exchange 
channel is discussed in Sec. 3.2. 



2.3 Application to Coulomb Blockade statistics 



A few experimental results on the statistics of the Coulomb blockade peak 
spacing in a quantum dot [7-9] suggest that the predictions of the constant 
interaction model, with single-particle levels taken from RMT and with a con- 
stant charging interaction U c only, fails to describe the fluctuations of Coulomb 
blockade peak spacing. The RMT+f/ c -model predicts a Wigner surmise peak 
distribution, and even-odd effects. But experimentally, the peak spacing dis- 
tribution is roughly Gaussian (with non Gaussian tails), it width is > A [9], 
and even-odd effects are not observed. Numerical studies for a few electrons, 
with mutual Coulomb repulsion, in a disordered medium [17] deviate from the 
predictions of the RMT+t/ c -model as well. 
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We will discuss now what the model (1) predicts for the Coulomb block- 
ade peak spacing distribution. By definition, the spacing between the iV'th 
Coulomb peak and the N — l'th Coulomb peak is given by 

AE = (E N+ i — E N ) — (E N — £jv-i) 

where E^ is the energy of the system with N electrons. Fluctuations in the 
single particle energy levels {s^}, and in the interaction parameters U c , J s and 
J c may lead to fluctuations in AE as iV changes. The latter, however, vanish 
when g — > oo. As we discussed in the previous section, the fluctuations of the 
single particle levels induce fluctuations in the ground state spin of the dot 
even when the interaction constants J s , U c and J c do not fluctuate with N. 
These, in turn, cause fluctuations in the peak spacing. 



Coulomb Peak spacing/ A 




Fig. 3. Theory for cumulative Peak spacing distribution in the presence of time 
reversal symmetry (/? = 1). We show the case without interaction in the Cooper 
channel (J c = 0) as it is suppressed (see Appendix C). For clarity we shift the 
distributions with exchange parameters A = — J s /A in intervals of half of the average 
level spacing, A. 

In Fig. 3 we plot the cumulative Coulomb peak spacing distribution for ensem- 
bles with a GOE symmetry. As discussed in Appendix C the interaction in the 
Cooper channel is suppressed (by a logarithmic factor). We therefore plot in 
Fig. 3 the cumulative peak spacing distribution for a GOE ensemble without 
interaction in the Cooper channel(J c = 0). Fig. 3 shows the distributions for 
exchange interaction strengths of A = — J s /A = 0,0.1, . . .0.8. We choose to 
plot the cumulative peak spacing distributions and not histograms of the peak 
spacing density distribution. Plotting the cumulative distribution allows us to 
represent delta functions in the spacing distribution and avoids the need to 
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work with arbitrarily binning intervals, as is needed for a histogram. 

Examining Fig. 3, one notices that there is a jump in the cumulative Coulomb 
peak spacing distribution at energy AE = 2AA, corresponding to a delta func- 
tion in the spacing distribution at that energy. This occurs because there is a 
finite probability that starting from a spin-singlet ground state, two successive 
electrons will enter with opposite spin into the same single-particle state, and 
the quantity 2AA is the exchange energy cost for adding the second electron, 
when J c = 0. However, at large values of A the probability that the ground 
state of the dot is a singlet is smaller, hence the hight of the jump decreases. 
There is an additional substructure of the distributions (e.g. a kink in the 
lower part of the distribution) that becomes smoother when A increases. To 
understand a few details of the Coulomb peak spacing distribution curve, we 
have to find the ground state energies of a disordered system with N — 1, 
N and N + 1 electrons. (We assume that the system has the same disorder 
realization for consecutive electron entries.) 

The ground state of a system with N electrons is also characterized by its 
spin S. Few examples of the energies of states \N, S) are summarized in Fig. 4. 
The Coulomb peak spacing is given by AE = £ , jv-i,s jv _i +En+i,s n+1 — 2£jv,s iV , 
where En,s — (N, S\ H, e g \N, S) is the energy of the state \N, S). [Notice that 
in the presence of interaction in the Cooper channel \N, S) is not necessarily 
an eigen state of H e g, defined in Eq. (1)]. Different values of SV-i, S N , S N+1 
give rise to different spin sequences. 

There are few possibilities for sequences of spin entries, (see Table 1). Sequence 
#1 describes a situation where initially there are (2m — 1) electrons in the 
system. The first 2(m — 1) single particle states are doubly occupied, and the 
last state (state m) is singly occupied. We denote this state by |l,y 2 ), see 
Fig. 4. Then an electron is added to state m so that it is doubly occupied, we 
denote this state by |2, 0). The third electron is added to state m + 1 so that it 
singly occupied to form the state |3, y^). The peak spacing of this sequence is: 

AEi = £1,1/2 + £3,1/2 — 2i?2,o 

= U C + % J s + 9U C + 5 + 3 / 4 J s + J c - 2(4C/ C + J c ) 
= 2U C + 3 / 2 J s -Jc + 5 = 5-u-J c = AT 1 . 

In the last equality we have used relation (3) for the contact interaction toy 
model. The peak spacings for other sequences, calculated in a similar way, are 
summarized in Table 1. 

Sequence #1 occurs only if £2,0 < £2,1 5 > J c — 2J S <^ A7\ > u. In a 
similar way one can check that sequence #2 occurs when AT2 > u. Thus, both 
processes #1 and #2 will lead to a step in the peak spacing distribution at u. 
For the cumulative peak spacing distribution, this will lead to a discontinuity 



11 



t J 



m+2 ' 
m+1 i 

m ■ 

m+2 ■ 
m+1 ■ 
m ■ 

m+2 ■ 
m+1 ■ 
m ■ 

m+2 
m+1 
m 

m+2 
m+1 
m 

m+2 
m+1 
m 

m+2 
m+1 
m 



|N,S> 



N,S 



■ 8 

■ § 
■0 

.5 
.8 

■ 



■8 

+1 



.8 
.8 

4-o 



.8 
.8 

^0 



0,0) 

lu/ 2 > 





IU 3 /u 



—s |2,0) 4U C+ J C 



—5 |2, 1> 4lU2J s + 5 



—s |3, 1 / 2 ) 9U c+ 3 / 4 J s + J c+ 5 



|4, 0) 16U c +4J c+ 25 



—5 |4, 1> 16lU2J s + J c+ 5+5 



Fig. 4. Possible spin configurations of the dot and their energies according to Eq. (1). 
We assume that all states below level m are doubly occupied and denote this many 
body state |0, 0). The single particle states are: e m = 0, e m+ i = 5, e m+ 2 = 8- 

in the slope of the curve, which may be seen in Fig. 3 at the points AE = XA = 
—J s . Processes #1 and #2 are the main contributions to the approximately 
linear portions of the curves in the range — J s = XA < AE < 2XA = —2J S , 
which are seen at small values of A. 



Sequence #3 will occur if _E 2 ,o < E 2> i and E i>0 < E^. This occurs if J C — 2J S < 
5 < 52 J s — 3J C and leads to a step function jump in the cumulative peak 
distribution at AE 3 . This is clearly seen in the theoretical curves of Fig. 3 and 
Fig. 5. The weight of the jump is bounded from above by J^ E p(S')dS' where 
p{5') is the Wigner distribution for consecutive levels at distance 5'. 
Sequence #5 requires that E 2i i < £2,0 and E^q < £4,1, which occurs when 



12 



1 


Spin Configurations 




AT- 


I 


1 1 1 / \ v 1 o n\ v 

|i, h) => u) 


1 9 

|3, 


1/ \ 
IV 




—7/ — 7 -I- A 





|1, V 2 > |2, 1) 


|3, 


1/ \ 

/2/ 


£U C 1 2<-'s i J c v 




3 


12, o) => |3, y 2 > 


4, 


0) 


2Ur. — 3 /o,L + ,L 


2u + Jr- 


4 


12, o) |3, y 2 > 


4, 


1) 


2f/ c + y 2 «/s + 5 


5 


5 


12, i) |3, y 2 > => 


4, 


0) 


2f/ c + 1 hJ s + 2J c + 5 


+ 2J c + 5 


6 


|2, 1) |3, y 2 > |4, 


1) 


2U C + 5 / 2 Js -Jc + S 


-2u -Jc + S 



Table 1 The symbol AEi denotes the spaces between the Coulomb peaks for spin 
sequences i. [ATj is AEi using the parameters of the toy model, Eq. (3).] Sequences 
that involve higher spins are also possible and are not included in this table. To find 
the actual contribution to the peak spacing we should include the probability that 
such a sequence occurs. 

5 < min{ J c — 2 J s , 5 — 3 J c + 2J S }. This leads to a low energy tail in the shape 
of the Wigner distribution curve, that extend all the way down to AE = 0, 
when J c = 0. 

Here we have discussed some of the simplest spin-entry sequences, including 
situations with total spin 0, 1/2, and 1, and we have seen how these sequences 
appear in the cumulative peak spacing distribution. Generalizations for more 
complex situations are straightforward but tedious. Different sequences will 
lead to other singularities and smooth curve-segments in the distribution of 
peak spacings. 

The actual numerical calculations (which includes also sequences with spin 
larger than 1) do not demand, however, a detailed analysis as above. We 
performed them in the following manner: we use 24 random levels around the 
center of the spectrum of a random 100 x 100 matrix for 1000 realizations. For 
each realization we find, using formula (4), the configuration of the groundstate 
and its energy for eight consecutive entries of electrons, this gives us 6 peak 
spacings for each realization, so that totally we plot a histogram of 6000 peaks. 

In Fig. 5 we plot the cumulative peak spacing distribution for magnetic field 
large enough so that the non interacting levels may be described by the GUE 
ensemble and the interaction in the Cooper channel is completely suppressed. 
(We assume that the magnetic field is small enough so we can neglect Zeeman 
splitting effects.) 

2. 3. 1 Comparison between Theory and Experiments 

This section compares our theory for Coulomb blockade peak spacing with 
the available experimental results. We will see below that the agreement be- 
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Fig. 5. Theory of cumulative peak spacing distribution in the absence of time rever- 
sal symmetry ((5 = 2). We assume that the external magnetic field is large enough 
to break time reversal symmetry, but small enough so that we can neglect Zee- 
man splitting. In other words, the interaction in the Cooper channel is completely 
suppressed and the level spacing distribution is described by a GUE ensemble. For 
clarity we shifted the distributions with exchange parameters A = — J s /A by inter- 
vals of half level spacing. 

tween theory and experiment is not very good, particularly for small peak 
spacings. However some features of the theory are found also in experiments. 
For example, a non-Gaussian tail at large values of peak spacing and a jump 
in the cumulative distribution, that we described in the preceding section, are 
present in few experiments. 

Fig. 6 depicts the cumulative peak spacing distribution of AlGaAs-dots. Dots 
"3-7" were studied in Ref. [9] and dot "S" in [10]. In all the experimental 
curves that we present here, no magnetic field is applied. We therefore assume 
in the theoretical analysis that time reversal symmetry is conserved. In each 
curve we normalize the peak spacing by the dot level spacing, which vary 
from dot to dot since their sizes are different. We have used here the average 
level spacings quoted in the experimental papers, there are, of course, some 
uncertainties in these values. After this normalization we would expect that the 
curves of dots 3-7 will be similar. This should happen because they have similar 
electron density and therefore similar r s and exchange interaction parameter 
(see Appendix B). In addition we would expect that dot "S" will have a 
different curve, because it has an electron density that is larger by a factor 
~ 3. However, as Fig 6 shows, the experimental results behave differently. This 
may be attributed to the intrinsic noise in the system (see Table 1 in [9]), but 
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Fig. 6. Experimental results (curves "3-7" [9] and "S" [10]) of the cumulative peak 
spacing distributions. To compare between the dots of different sizes, we subtract 
from each distribution the dot average peak spacing and scaled it to the average 
level spacing, A. The dots here have time reversal symmetry. 



we still do not understand completely the origin for this behavior. The data of 




Fig. 7. Comparison between theory and experiment [10]. To fit best the upper 
tail of the experimental results in the absence of magnetic field, we choose in 
the theory (with GOE symmetry but without interaction in the Cooper channel) 
-J s /A = A = 0.4. 

Refs. [7] and [8] show a significantly wider distribution of level spacings than 
found in Ref. [9], despite the apparent similarity of the systems studied by 
these groups. So far, there has been no clear explanation for this discrepancy. 



We nevertheless plot in Fig. 7 our theory and the experimental curve of [10] 
(without magnetic field). The overall experimental curve fits well to a Gaus- 
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sian. However, the details of the upper (right) tail of the cumulative distri- 
bution, i.e., the jump and the non Gaussian tails fits better to the theory 
of the spin fluctuations in the Ground state. To fit best the upper tail of 
the experimental results in the absence of magnetic field, we choose in the 
theory (with GOE symmetry but without interaction in the Cooper channel) 
— J s /A = A = 0.4. Notice that using the static RPA estimate of A (in Ap- 
pendix B), with the experimental value r s ~ 0.72 [10] we find, as expected, a 
value that is somewhat smaller than 0.4. 

We expect that several effects, (not included in the theory) may be important 
for current experiments. The lower part of the peak spacing distribution, that 
is built from single particle levels that are, by chance, very close to each other, 
is especially sensitive to these effects. Indeed, this part appears to be far from 
the theoretical curves. Among them are: 

(1) Non universal effects of finite g that cause fluctuations in the interaction 
parameter [12]. For ballistic two-dimensional dots g = \ / 2rmA, where 
A is the dot area and n is the density of the electrons. In the dots of 
Ref. [7-9] r s ~ 1 - 3 and g ~ 50 - 150 (in these of Ref. [10] r s ~ 0.72 and 
g ~ 35). Thus, when we compare theory to experiment, fluctuations in 
the interaction constants cannot always be ignored. Ullmo and Baranger 
present in Ref. [12] a detailed study of the effect of fluctuations of the 
interaction parameters (see also the discussion here in page 6). They 
find indeed that when these fluctuations are included the results "are 
significantly more like the experimental results than the simple constant 
interaction model" . 

(2) We assume that the temperature, and the single particle levels width 
(due to tunneling to the leads), is smaller than the mean level spacing 
and therefore we ignore their effects. This assumption is not valid for the 
lower part of the distribution as it is built from levels whose distance from 
neighboring levels might be much smaller than the average level spacing. 
The importance of the temperature was considered very recently by Usaj 
and Baranger [22] . They find that temperature effects are significant even 
at T ~ 0.1A. 

(3) There is experimental noise due to charge motions during the measure- 
ments time. This effect [9] might be the dominant contribution to the 
smearing of the distribution in the experimental curve. 



3 Spin-Orbit Effects 

Spin-orbit coupling can have major effects on the ground states or the low- 
energy transport properties of a mesoscopic system. In many metallic nanopar- 
ticles, spin-orbit effects arise from randomly placed heavy-ion impurities, which 



16 



can simultaneously scatter electrons and flip their spin, subject to the con- 
straints imposed by the requirement of time-reversal invariance in the absence 
of an applied magnetic field. In other cases, one is concerned with metal parti- 
cles where the spin-orbit effects are already significant in the band-structure of 
the ideal host crystal, so that the "spin" variable in the Bloch states actually 
represents a mixture of spin and orbital degrees of freedom at the microscopic 
level. In this case spin-flip scattering with the requisite spin-orbit symmetry 
can occur whenever there is scattering: from defects, from impurities, or from 
the boundaries of the sample. Spin-orbit scattering in the above cases can 
generally be characterized by a spin-orbit scattering rate, and the importance 
of spin-orbit effects is determined by the ratio of this rate to other frequencies 
characteristic of the mesoscopic system. Effects of spin-orbit scattering on the 
groundstate spin-structure and on the energy splitting in an applied magnetic 
field will be discussed in Subsection 3.1. The effects of spin-orbit coupling 
on the spacing of groundstate energies, in the presence of electron-electron 
interactions, will be discussed in Subsection 3.2. 

A peculiar situation can arise in two-dimensional electron systems in materials 
such as GaAs. Here the dominant spin-orbit effects arise from terms in the 
effective Hamiltonian in which there is a coupling to the electron spin linear in 
the electron velocity. (These terms arise from the asymmetry of the potential 
well confining the electrons to two dimensions and from the lack of inversion 
symmetry in the GaAs crystal structure.) The special form of this coupling 
leads to a large suppression of spin-orbit effects when the 2D electron system 
is confined in a small quantum dot. However, effects of spin-orbit coupling are 
again enhanced in the presence of a strong magnetic field parallel to the plane 
of the sample, so that they must be taken into account in such properties as 
the level-spacings of a closed dot or the statistics of conductance fluctuations 
in a dot coupled to leads through one or more open channels. These effects 
will be discussed in Subsection 3.3 below. 



3.1 Effective g-tensor of a metal particle with spin-orbit scattering. 

According to Kramers theorem, a metal particle with an odd number of elec- 
trons with no special symmetry, in zero magnetic field, must have a degener- 
ate groundstate manifold, with pairs of states related to each other by time- 
reversal symmetry. In the absence of spin-orbit coupling, the total spin S is a 
good quantum number, and the groundstate manifold is just that expected for 
half- integer S. As we have seen in Sec. 2, if the electron-electron interaction is 
weak, we will essentially always find S — 1/2 for odd iV and the ground state 
will be just two-fold degenerate. For stronger electron-electron interactions, 
however, there will be some probability of finding S = 3/2 or larger, so that 
four-fold or higher degeneracies are also possible. When spin-orbit interactions 
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are turned on, the higher degeneracies will be broken into a set of doublets, 
so that the ground state will again be two-fold degenerate. 

If we now apply a magnetic field B to the system, the degenerate ground state 
will be split. For sufficiently small B, one of the states will move up in energy 
by an amount Se which is linear in B, while the other state will move down 
by the same amount. These shifts may be measured, at least in principle, 
by electron-tunneling spectroscopy experiments in an applied magnetic field. 
We discuss here the statistical properties of the distribution of energy shifts 
expected under various circumstances. We concentrate on the situation where 
the electron-electron interaction is weak, so that the many-body ground state 
is well described by the picture of weakly interacting quasiparticles, as effects 
of electron-electron interactions will be discussed in the next subsection. 

Quite generally, we may write the linear splitting of a Kramers doublet in the 
form 

Se= \fi B /2\(B ■K-B) 1 ' 2 (5) 



where //# < is the electron Bohr magneton and K is a real, positive-definite 
symmetric 3x3 tensor. In the absence of spin-orbit coupling, K is isotropic, 
with Kij = 45ij. When spin-orbit coupling is present, we find that K varies 
from level to level, and is in general anisotropic. We write the three eigenvalues 
of K as gl, (k = 1,2,3), with \gi\ < \g 2 \ < \gs\, and refer to the g^s as the 
three principal g-factors for the level. Although the energy-splittings in a static 
magnetic field only define the absolute values of the gk, by considering the 
response to a time- varying magnetic field (e.g., a spin resonance experiment) 
one can also give an unambiguous meaning to the sign of the product of the 
three ^-factors. Since the sign of an individual gk has no physical meaning, we 
adopt the convention that g 3 and g 2 are always positive, but g 1 can be positive 
or negative, depending on the specific system considered. 

For the case of weakly interacting electrons, which we consider here, the ground 
state for 2N + 1 electrons consists of 2N electrons in filled Kramers doublets, 
plus one electron in a doublet which is singly occupied. The filled doublets 
give no contribution to the linear energy shift because in each case one state 
moves up and the other moves down by the same amount. Thus, the (/-factors 
are determined by the properties of the singly-occupied state. 

In the presence of spin-orbit coupling there are two contributions which can 
shift the g- values from the bare value g — 2. If we take into account only the 
interaction of B with the electron spin, then spin-orbit coupling will always 
reduce the ^-values. For example, if the magnetic field is applied in the z- 
direction, the state which is shifted down in energy will be the particular linear 
combination of the two degenerate states which has the maximum expectation 



18 



value of — S z . This expectation value is < 1/2, so the spin-contribution to the 
g-factor will generically be reduced by spin-orbit coupling. 

On the other hand, there is also an orbital contribution to the linear Zeeman 
effect, when spin-orbit coupling is present . (In the absence of spin-orbit cou- 
pling, the orbital states in an irregular dot will be generically non-degenerate 
and time-reversal- invariant, so they cannot acquire a linear energy shift in a 
weak magnetic field.) Both the orbital and spin contributions were consid- 
ered by Matveev et al. [23], who discuss the expectation value and probability 
distribution of Se 2 for the magnetic field in an arbitrary fixed direction. 

By contrast Brouwer et al. [24] considered the joint probability distribution 
of the three g-values for a single level, so they could examine the anisotropy 
as well as the magnitude of the g-tensor. Their analysis concentrated on the 
case where orbital effects can be ignored, so that the mean-square g-factors 
are monotonically reduced with increasing spin-orbit coupling. The strength 
of the spin-orbit coupling in this case is determined by a parameter 

where A is the mean separation between one-electron energy levels. The mean 
spin-orbit scattering time r so is defined so that if we prepare a state with 
spin up, the probability to find it in the same spin direction after time t is 
~ e -*/ r s°. When A so 3> 1, one finds that the g-factors are greatly reduced from 
their bare values, and one can obtain an analytic form for the joint probability 
distribution: 

P(gi,g2,g 3 )<xl[\g*-g*\l[e~^^ 2 ) (7) 

i<j i 



For intermediate values of the coupling parameter A so , one can perform nu- 
merical simulations to study the distribution, using random-matrix theory. In 
Fig. 8 we show the A so -dependence of the mean values of gf, as well as the 
values of gu for a particular realization of the random matrices. 

Very recently, Petta and Ralph [Ref. [25]] have measured effective g- values for 
a number of levels in each of several different nanoparticles, of Cu, Ag and 
Au, with diameters in the range 3 - 5nm. They did not vary the direction 
of the applied magnetic field, so they could not study the anisotropy of the 
g-tensor. However, the statistical distributions of the g-factors (normalized to 
their means) for different levels in a given particle were found to be in good 
agreement with the theories of Refs. [23] and [24]. For the mean g-factor, an 
agreement with Refs. [23] and [24] is found if the spin contribution is taken into 
account only; the mean g-values observed in the Au particles (ranging from 
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^so 

Fig. 8. Average of the squares of principal ^-factors versus spin-orbit scattering 
strength A so , obtained from numerical simulation of a random matrix model. Inset: 
<7i, g2, and g% for a specific realization. We have included the sign of g\ . 

0.12 to 0.45) were significantly smaller than what one might expect from the 
orbital contribution, according to the theory of Ref . [23] , unless one assumes 
a very short mean-free-path for the electrons. (Using the formulas in Ref. [23], 
one would need a mean-free-path of order 0.1 wn to get g-values this small.) 
Small g-values (below 0.5) for Au nanoparticles were also observed previously 
by Davidovic and Tinkham [4]. 

3.2 Effects of spin-orbit coupling on interaction- corrections to groundstate 
configurations and energy spacings. 

So far, we have ignored the effects of electron-electron interactions. This should 
generally be valid if the exchange interaction is small compared to the thresh- 
old for the Stoner instability, so that the probability of finding S > 1/2 is small 
in the absence of spin-orbit coupling. When the spin-orbit coupling parameter 
A so is large, the effective exchange interaction between two electrons in states 
close to the Fermi-energy of the particle will be even further reduced, as the 
mean spin in any state becomes small compared to 1/2, and the local spin ori- 
entations have different spatial distributions for one-electron states belonging 
to different Kramers doublets. 

In the limit of very large spin-orbit coupling, where the mean spin tends to 
zero, the exchange interaction should also tend to zero. This means that the 
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parameter J s in the effective Hamiltonian (1) of Sec. 2.1 should be set to zero. 
This is consistent with the fact that spin is no longer a good quantum number 
of the system, and the term proportional to J s is no longer invariant under 
the set of allowed unitary transformations of the random matrices. 

A consequence of this analysis is that if a spin-orbit scatterers are added 
to a system with fixed electron-electron interaction (fixed r s ) the probability 
distribution for the separation of successive groundstate energies, measured 
by the Coulomb-blockade peak separations, should approach that of a non- 
interacting electron system in the symplectic ensemble. This means that there 
should be a bimodal distribution with an even-odd alternation. The chemical 
potential to add a second electron to a Kramers doublet is the same as the 
energy to add the first electron, after the coulomb blockade energy [U c in 
Eq. (1)] is subtracted, whereas the chemical potential for the next electron 
will be larger by an amount approximately given by the mean level spacing A. 

3.3 Spin-orbit effects in a GaAs quantum dot in a parallel magnetic field. 

The most important spin-orbit terms in the effective Hamiltonian for a 2D 
electron gas (2DEG) in a GaAs heterostructure or quantum well may be writ- 
ten in the form 

Ti so = 7iv x cr y - ^ 2 v y o- x (8) 

where v is the electron velocity operator. We have assumed that the 2DEG 
is grown on a [001] GaAs plane, and we have chosen the x and y axes to lie 
in the [110] and [110] directions. For an open 2DEG this leads to a spin-orbit 
scattering rate of order 7 2 -D, where D is the diffusion constant and 7 is the 
geometric mean of the two coupling constants in Eq. (8). For a confined dot 
of radius R, in zero magnetic field, however, the effects of spin-orbit coupling 
are suppressed if the typical angle of spin precession for an electron crossing 
the dot, given by 9 = 7/?, is small compared to unity. One finds in this case 
that the matrix elements of Ti so are greatly reduced for energy states whose 
energy separation is small. 

Halperin et al. [26] have argued that the effects of spin-orbit coupling can be 
enhanced, however, in the presence of an applied magnetic field in the plane. 
The enhancement is maximum when the Zeeman energy becomes comparable 
to the Thouless energy (i.e., the inverse of the transit time for an electron in 
the dot), in which case there is an effective spin- mixing rate comparable to 
the spin-orbit scattering rate for an open system with an electron mean free- 
path equal to the mean free path in the dot. For a closed dot, the spin- mixing 
would be manifest in the repulsion of energy levels for different spin, and the 
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appearance of anti-crossings of the levels as when the Zeeman field is varied. 

Motivated by experiments of Folk et al. [27], Halperin et al. [26] considered 
the "universal conductance fluctuations" of a dot connected to a pair of leads 
with one or more channels open in each lead. They considered explicitly the 
case where there is a weak magnetic field perpendicular to the dot, so that 
time reversal symmetry is broken, and the system is in the class of the unitary 
ensemble, even in the absence of spin-orbit coupling. It was shown that effects 
of spin-orbit coupling in large Zeeman field could then lead to a factor of two 
reduction in the variance in the conductance, which is in addition to the factor 
of two reduction caused by breaking of the spin degeneracy. Calculations of 
the cross-over, as a function of the in-plane magnetic field, were in at least 
qualitative agreement with the experimental observations. 

Very recently, Aleiner and Fal'ko[28] have considered the case without a per- 
pendicular magnetic field, so that the system without spin-orbit coupling 
would be in the orthogonal ensemble. They have shown that the application 
of a parallel magnetic field in this case turns on a spin-orbit perturbation with 
a special symmetry, so that the system retains an effective time-reversal sym- 
metry even in the presence of the large Zeeman field. The spin-orbit coupling 
leads to a reduction in the size of conductance fluctuations, but not as much 
as one would obtain if the time-reversal symmetry was also broken. The spin- 
orbit coupling also leads to a reduction in the "weak localization" correction 
to the average conductance, but does not lead to complete suppression as one 
would find for a broken time reversal symmetry. (However, as noted by Meyer 
et al. [29] and by Fal'ko and Jungwirth[30], for an asymmetric quantum well 
of finite thickness, application of a strong magnetic field parallel to the sample 
can lead to broken time-reversal symmetry due to orbital effects, even in the 
absence of spin-orbit coupling.) 



4 Origin of multiplets in the differential conductance 

In a recent experiment Davidovic and Tinkham[4] studied tunneling into indi- 
vidual Au nanoparticles of estimated diameters 25nm, at dilution refrigerator 
temperatures. The differential conductance dl/dV, as a function of the source- 
drain voltage V, indicate resonant tunneling via discrete energy levels of the 
particle. Unlike previously studied normal metal particles of Au and Al, in 
these samples they find that the lowest energy tunneling resonances are split 
into clusters of 2-10 sub-resonances. The distance between resonances within 
one cluster is much smaller than the mean level spacing of the Au grain. 

This situation is illustrated schematically in Fig. 9. The differential conduc- 
tance dl/dV shows resonances, where each resonance in dl/dV is actually 
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a multiplet, the splitting between the peaks of the multiplets being a factor 
~ 30 smaller than distance between the resonances (which is of the order 
single-particle level spacing in the grain). 

In this section we outline two-different mechanisms which can lead to a fine 
structure of the first conductance peak. We first show how such a fine structure 
can occur if the ground state has a finite spin with small energy splittings 
between states of different magnetic quantum number. In this model it is 
necessary to have a relatively large total spin in order to split the conductance 
peak into many sub-peaks. This mechanism would also be suppressed by large 
spin-orbit coupling. 

In the second mechanism, following Agam et al. [31], we show how such a fine 
structure can arise from nonequilibrium processes induced by the large bias 
voltage V used in the experiment. This mechanism seems to us to be the more 
likely one for explaining the observations of Ref. [4]. We also indicate how, 
experimentally, one might distinguish between the two proposed explanations. 




V 

Fig. 9. Schematic illustration of the experimental results of Ref. [4]: The peaks in 
the differential conductance are split. The distance between the multiplets is of the 
order of the single particle level spacing A; the distance between peaks in the same 
multiplet is much smaller. 



4-1 Multiplets from an almost degenerate groundstate 

In general, a peak in the differential conductance as a function of the bias 
voltage V may occur if an additional channel for tunneling onto or from the 
metal grain is opened at that V. The relation between V at the peaks and 
the ground state energies is complicated; it depends on the capacitive division 
between the left and right contact and on the conductances of the two tun- 
neling contacts. A detailed account of the possible scenarios can be found in 
the review by von Delft and Ralph [32]. 

Here, we make the simplifying assumption that the left point contact has 



23 



the bigger resistance and the smaller capacitance, so that the electrostatic 
potential of the dot equals that of the right reservoir, and the contact to the 
left reservoir can be seen as the "bottleneck" for current flow. Then, if the 
grain has N electrons at zero bias, a conductance peak occurs when 

eV = E N+1 — E N , 

i.e., when the bias voltage V is precisely equal to the difference of the energies 
of any two many-body states of the grain with iV and N + l electrons (for 
V > 0), provided the initial iV-particle state is populated at the corresponding 
bias. Below we focus on the first peak in the differential conductance and 
discuss when and how a fine structure of that first peak can arise. We assume 
that the temperature is small compared to any splittings in the energy levels, 
and we assume (for the moment) that there is no spin-orbit coupling. 




Fig. 10. Schematic drawing of the tunneling process considered here. The left point 
contact has the smaller capacitance and the smaller conductance. When the bias 
voltage is increased, peaks in the differential conductance occur, whenever a new 
channel for tunneling onto or from the grain is opened. Compare the bias voltages 
in (a) and (b). 

If the N and iV + 1-particle ground states have perfect degeneracies, the dif- 
ference E N — E N+1 can only take a single value, and a single peak will be 
observed, no matter what the spin of the ground state is. Hence to observe a 
fine structure, the degeneracy of the ground state has to be lifted. This can 
be done by application of a uniform magnetic field, as is illustrated in Fig. 11 
for Sn = or Sn = 1 and Sn+i = 1/2. In the case where Sn+i = Sn + 1/2, 
the difference E^ + i — E^ between the energies of the many-body states for 
iV and N+l particles can take two values, 



E 



N+l 



-E 



N 



E% +1 -E N ±{l/2) 9l i B B, 



where E% and E° N+l are the N and (N + l)-particle energies in the absence 
of the magnetic field. The differential conductance shows a double peak at 
voltages 



eV+ 



E 



N+l 



E N ± (l/2) gf i B B, 



as is seen in Fig. 11a. On the other hand, only a single peak at bias voltage 
V + = E N+ i — E N + (l/2)g/i B B is found if SV+i = Sn — 1/2. Although the 
bias voltage V_ corresponds as well to a transition energy between many-body 
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states with N and N + 1 particles for Sn > Sjv+i, no peak in the differential 
conductance is found at that bias voltage, because the initial state of that 
transition is an excited state, which is not populated at V = V_. Population 
of an excited iV-particle state is only possible at higher bias voltages V >V + 
via inelastic processes that use the N + 1-particle state as an intermediate 
step. (A small nonequilibrium population of the excited iV-particle state, and 
hence a small peak at V = VI, may, however, occur as a result of inelastic 
cotunneling, as is explained in Ref. [33].) If the difference in the total spin 
quantum numbers for N and iV + 1 is greater than 1/2, then there can be 
no conduction peak at all in the absence of spin-orbit coupling or inelastic 
cotunneling processes. 




(a) S N =0 S N+1 =l/2 S N =1 S N+i = 1/2 (b) 

Fig. 11. Possible transitions between Zeeman split states with N and N+1 particles, 
for Sn = 0, Sn+i = 1/2 (a), and Sn = 1 and Sjv+i = 1/2 (b). Note that the 
transitions starting out of the excited states of the triplet in (b), denoted by the 
dashed arrows, do not give rise to peaks in the differential conductance, (assuming 
that equilibrium is reached between successive tunneling events) , because the excited 
states are not populated at eV = En+i — En — fJ-BgB/2. 

The situation changes in the presence of weak static magnetic impurities. 
"Weak" here means that they can be seen as a small perturbation on top of 
the picture sketched in the previous sections. "Static" means that the impu- 
rity ion has a large intrinsic angular-momentum and large crystal-anisotropy, 
so that we can neglect the matrix elements for transitions between different 
impurity spin-states. The impurity spins could be in the grain itself or could 
be located close to the grain in the surrounding insulator. Then if the many- 
body state of the grain has non-zero spin, the spin degeneracy will be lifted 
by the coupling to the impurity spin, which can give rise to a splitting of the 
lowest conductance peak even in the absence of an applied magnetic field. A 
significant difference between this case and the splitting due to an external 
field is that the effective coupling now depends on the microscopic details of 
the electron wavefunctions close to the impurity, and the level splitting will 
generally be different for the N and N+1 electron states. As we shall see, 
this makes it possible for the lowest conductance resonance to split into more 
than two sub-peaks. 

We first consider the case of a single impurity spin. According to the Wigner- 
Eckart theorem, if the coupling to the impurity spin is weak, an electronic 
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many-body state with total spin S will be split into (25* + 1) equally-spaced 
levels, characterized by the quantum number of the magnetic moment in the 
direction parallel to that of the frozen impurity spin. The size of the splitting 
depends on the concentration and microscopic details of the impurities. Since 
a peak in dl/dV can occur whenever eV = E^+i — En, many close peaks 
appear when the degeneracy of the ground state is lifted. The total number of 
possible transitions is (2Sn + 1)(2Sjv+i + 1), since now En and En+i can take 
2Sn + 1 and 2Sn+i + 1 values, respectively. However, for the same reasons as 
discussed above, not all possible transitions give rise to peaks in the differential 
conductance: Only transitions at energy differences AE = E N+ i — E N where 
the initial iV-particle state is already populated at a bias voltage eV < AE 
are reflected as peaks in the differential conductance, and the spin component 
parallel to the frozen impurity spin can only change by ±1/2. Some examples 
are shown in Fig. 12 for S N = 1/2, S N+ i = 1 and S N = 1, S N +i = 1/2. In 
the figure, the transitions that correspond to true peaks in dl/dV are shown 
as solid arrows, the other ones are shown with dashed arrows. In practice, 
since eV is typically much bigger than the fine structure of the N and N + 1- 
particle levels, all transitions appearing at energy differences AE bigger than 
the difference el^h = E N+1 — E 9 N between the ground state energies for N 
and N + 1 particles will show up as true peaks at eV = AE , while no peaks 
appear for eV < eV t h (V t h is the threshold voltage for current flow). 




(C) N N+l N N+l (d) 

Fig. 12. Possible transitions between states with SV = 1/2, SV+i = 1 (top) and 
Sn = 1, Sn+i = 1/2 (bottom), when the spin degeneracy is broken by the presence 
of (several) static magnetic impurities. The transitions indicated by solid arrows give 
rise to peaks in the differential conductance; the transitions indicated by dashed 
arrows do not, because the state they are starting from is not populated at the 
corresponding bias voltage. How many peaks are visible depends on the actual 
splitting of the energies. 
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If there are several frozen impurity spins coupling to the electrons, we can 
again use the Wigner-Eckart theorem, in the case of weak coupling, to show 
that a groundstate with spin S is split into (2S + 1) equally spaced levels 
classified by the spin component along some direction which is a weighted 
vector sum of the several frozen impurity spins. The weights in this sum will 
be different in the states with N particles and with N + 1 particles, so that 
the quantization axes will generally be different in the two states, as well as 
the level splittings. This will lift the selection rule that the quantum number 
can only change by ±1/2 on the addition of a single electron, so the number 
of lines in the multiplet may increase accordingly. 

The model of coupling to one or more frozen spins can also be generalized to 
the case of dynamical spins. For example, in the case of coupling to a single 
dynamical localized spin Si in the insulating material close to the grain, the 
localized spin and the spin of the electrons in the metal grain S together will 
form states with total angular momentum ranging from \S — Si\ to |S'+S'j|, and 
will split up in the corresponding multiplet. In this case the number of possible 
transitions depends on the detailed selection rules governing transitions of the 
localized spin. 

Finally, we consider the situation where spin-orbit coupling is present. In the 
case of weak spin-orbit coupling, a groundstate with spin greater than 1/2 
can be split into several different energy levels even in the absence of an 
applied magnetic field and in the absence of magnetic impurities (although 
the splitting by spin-orbit coupling only arises in second order perturbation 
theory, whereas the splitting caused by magnetic impurities already appears 
in first order perturbation theory). In general, the various states will be split 
by different amounts, and so multiple subpeaks can be observed, for large S, 
just as we found for the case with a frozen magnetic impurity. In the present 
case, however, states with odd N remain twofold degenerate by Kramers' 
degeneracy, which reduces the number of possible transitions roughly by a 
factor of two, compared to the case of static magnetic impurities. 

If the spin-orbit coupling is too strong, however, spin-orbit splittings will be- 
come comparable to or larger than the single-particle level spacings. In this 
case, it is no longer possible for splittings between spin states to give rise to 
fine structure of the conductance peak on a scale small compared to A. (Also, 
as we have seen previously in Subsection 3.2, exchange-splittings tend to be 
reduced in this case, so that the one electron picture should be valid for the 
ground states.) The gold particles studied in Ref. [4] appear to be in the strong 
spin-orbit coupling regime. 
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4-2 Multiplets from nonequilibrium processes 

A second mechanism to observe multiple peak structures in the differential 
conductance is via nonequilibrium population of highly excited states of the 
metal grain, as was first proposed by Agam et al. [31]. This mechanism does 
not need a degeneracy, or near-degeneracy, of the ground state. The idea of 
Ref. [31] is as follows: Since the bias voltage is typically much larger than the 
spacing A between single particle levels, after an electron has tunnelled on 
and off the grain, the grain may be left in an excited iV-particle state, with 
an occupation of the single-particle levels that differs from the ground state, 
see Fig. 13. The fact that there is a different occupation of the single-particle 
levels will slightly shift the addition energies E N+ \ — E N , thus giving rise to 
peaks in dl/dV at different values of the bias voltage V. In Ref. [31] spinless 
particles were considered. In that case, nonequilibrium processes cause a fine 
structure of the second and higher resonances, but not the first one [31]. For 
spin 1/2 particles and even N, the scenario of Ref. [31] can also lead to a fine 
structure for the first resonance, as we will now describe. 

We denote the highest occupied (self-consistent) single-particle level in the TV- 
particle ground state by En/2 and assume that the iV-particle ground state has 
zero total spin. When the bias voltage exceeds the threshold eV t h = E g N+1 —E 9 N) 
current flow can leave the grain in an excited iV-electron state, when, after 
an electron has tunnelled into the level £jv/2+i, another electron tunnels out 
of a lower-lying level e U) see Fig. 13. Note that the excited state can have 
a total spin S — or S — 1. Since the grain is now in an iV-particle state 
that is different from the ground state (compare Figs. 13a and c), the energy 
cost Ejy + i — En for addition of an electron, and hence the position of a peak 
in dl/dV, is, in general, different from E 9 N+1 — E g N . A priori, this difference 
can have three contributions: (1) An electron can tunnel into different single- 
particle levels than in the ground state. (2) The transition energy E N+ i — E N 
depends on the spin S of the states involved, which can be different from the 
ground state spin. (3) All transition energies depend uniquely (but weakly) 
on the populations of the initial and final states through mesoscopic fluctu- 
ations of the interaction contribution to the energy [31]. The characteristic 
energy scales for the first two of these contributions are A and J s , while the 
third is of order A/^/g, as we shall see below, where g = 2-kE^/A is the 
dimensionless conductance of the grain. For an even N, however, there exist 
excited iV-particle states for which the level £jv/2+i is only singly occupied 
and S = 0, so that the first two contributions vanish. Then only the contribu- 
tion from mesoscopic fluctuations remains. If g is large, the energy scale for 
the fluctuations is small, and one finds multiple conductance peaks close to 
eV = E 9 N+1 — E 9 N , where the total width of the multiplet is of order A/y/g. 

Although all nonequilibrium configurations have their own characteristic tran- 
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Fig. 13. A nonequilibrium configuration can be obtained from the ground state (a) 
if an electron tunnels into the first empty level (b), and another electron tunnels off 
the grain from a different, lower lying, level (c). The energy required for tunneling 
an electron into the highest level in the configurations (a) and (c) may be different, 
which explains that more than one peak can be seen in the differential conductance. 

sition energy E^+i — En, not all of them need to correspond to a peak in the 
conductance; no peak is observed if the corresponding voltage is below the 
threshold Vth, which was needed to populate the corresponding excited state. 
We estimate that the number of sub-peaks in the first peak mutiplet, iV first _p Cak 
due to the nonequlibrium effect is of order eV/2A, which is roughly the ratio 
of the Coulomb blockade energy to the single-particle level spacing. 

To understand the reasons for this estimate let us first consider the case where 
we can neglect the exchange coupling J s in the effective Hamiltonian (1). (This 
is correct in the case where spin-orbit coupling is strong). We also neglect 
the Cooper-channel interaction J c , for the reasons discussed in Appendix C. 
However, we take into account fluctuations in the size of U c between different 
pairs of levels. In this case, the many-body states can be labelled by the 
occupancies of the single-particle levels, and the ground states with energy E 9 N 
and energy E g N+l are described by Fig. 13a and b respectively. Let us assume 
that V is close to, but slightly above, the threshold voltage Vth = E 9 N+1 — 
E g N . We further assume that when an electron enters or leaves the system 
it does not excite other electrons by multi-electron processes; in particular 
we ignore Auger-like processes. Similarly we assume that the only important 
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contributions to the conductance are via real, energy conserving, transitions. 

Under these assumptions, one finds that the excited states that contribute 
to iVfirst-peak have precisely one hole below En/2- The energy distance between 
these single-hole states is roughly A. Since only excited states within an energy 
eV from the ground state can be populated, the number of possible hole 
states is ~ eV/A. Of these, roughly half will lead to positive energy shifts, 
which are necessary for contributions to the multiplet structure, so we find 
IVfirst-peak ~ eV/2A. The width of the peak is proportional to the size of the 
fluctuations in U c that we have discussed in Subsection 2.1.1; in case of long 
range Coulomb interaction it is ~ A/^fg. 

Many-body states with two or more holes below e N /2 cannot contribute to the 
conductance for voltages close to V th , because if two holes are present when 
there are N electrons on the particle, the level e N/2+1 will necessarily be doubly 
occupied. Then, the next electron would have to enter through the level e N/2+2, 
which would require an additional energy, of order A, that is not available for 
V slightly above Vth- If, by chance, the distance of level e N/2+2 from level 
e N/2+1 is smaller than J s then transition through nonequilibrium states with 
two holes may occur. This situation will increase iV first _ peak significantly, it 
should be proportional now to ~ (U/A) 2 /2. 

Let us now consider the case where the exchange parameter J s is not negligible. 
Then, the N electron states having one electron in the level e N /2+i and one 
hole below the Fermi energy will have different energies depending on whether 
they have S — or S — 1. For S = 1, the energy is reduced by | J s \, so that 
the energy to add the next electron to the level e N/2+1 is increased by |J S |, 
and the peaks arising from the triplet configuration will be shifted up by this 
amount relative to the singlet contributions. If this shift is comparable to the 
level spacing A, then only the singlet peaks will appear close to the threshold, 
and the number of peaks in the lowest multiplet will be the same as before, 
-^first-peak ~ eV/2A. If \J S \ is sufficiently small, but not negligible, however, 
the singlet and triplet peaks may appear to form a single multiplet, with twice 
as many peaks as before. 



4-3 A comparison between the mechanisms of Subsections 4-1 and 4-2 

The main difference between the two explanations offered here is that the for- 
mer entails many different transitions between states very close to the N and 
N + 1-particle ground states, whereas the latter makes use of, in principle, 
the same transition between different pairs of states that are highly excited 
above the ground state. Thus, one may distinguish between the two scenarios, 
when the electrostatic potential of the grain can be changed by the voltage 
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Vg on a nearby gate: Fine tuning of V g affects the threshold bias voltage V t h, 
and hence the number of possible nonequilibrium configurations. Hence, if the 
fine structure of the first resonance is due to nonequilibrium processes, the 
peaks will disappear one by one when V g is tuned to the charge degeneracy 
point. This is illustrated in Fig. 14. In this figure, we have simulated the dif- 
ferential conductance from the rate equations of Ref. [5,34], for the case where 
all relaxation of excited states inside the grain occurs due to the coupling to 
the leads. The four panels show the differential conductance for four differ- 
ent values of the gate voltage, where the multiplet consists of 5, 4, 3, and 2 
peaks. The closer the gate voltage is to a charge degeneracy point, the fewer 
nonequilibrium peaks can be observed. On the other hand, if the fine structure 
is due to a degenerate ground state, no highly excited states are involved, and 
the number of peaks in the multiplet is insensitive to V g . Of course, both ex- 
planations (nonequilibrium processes and an almost degenerate ground state) 
can apply at the same time. In that case, multiple peaks will disappear at the 
same time, when V g is tuned to the charge degeneracy point. An alternative 
way to distinguish between the two scenarios is if the parity of the number 
of electrons iV can be changed by a gate voltage. Nonequilibrium processes 
cannot explain a fine structure of the first resonance if N is odd. 



5 Conclusions 

In this article, we have reviewed several effects related to the electrons' spin 
in the presence of spin-orbit coupling and/or electron-electron interaction, in 
small quantum dots with relatively large number of electrons. 

At low temperatures, in the absence of spin-orbit coupling, and in the absence 
of an applied Zeeman field, our system can be described by an effective Hamil- 
tonian of the form (1). In the limit where the number of electrons in a chaotic 
dot is large, the effective Hamiltonian contains three interaction parameters, 
in the direct, exchange and Cooper channels. In Appendices B and C we esti- 
mate these parameters for realistic systems. We have used use them in Sec. 2 
to calculate (i) the probability for the dot to have a non-zero total spin in its 
ground state and (ii) the distribution of the Coulomb blockade peak spacings. 
The theory for the latter describes many features of the experimental obser- 
vations, and is qualitatively much better than what one would have obtained 
if one ignored the exchange interaction. However, the simple effective model 
does not do well in describing the low energy tail of the distribution, and it 
does not account for the large differences in the data obtained by different 
experimental groups. 

In Sec. 3, we reviewed some recent theoretical studies on the effect of spin-orbit 
coupling in a quantum dot or metal particle. In Subsection 3.1 we discussed 
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Fig. 14. When the voltage V g of a nearby gate is varied, the number of excited 
TV-particle states that can be populated by current flow is also changed. The four 
panels show how the peaks in the differential conductance disappear one by one 
when V g is tuned in the direction of the charge-degeneracy point (at which current 
flow happens at zero bias voltage). As seen in the figure, the (minimal) bias voltage, 
V, for nonzero current flow decreases as V g approaches the charge-degeneracy point. 
In the four panels shown, the number of peaks decreases from five (upper left) to 
two (lower right panel). The single-particle levels are for the N + 1-particle system; 
the arrows indicate from which levels electrons can escape to the right reservoir. The 
dl/dV graphs were calculated using the rate equations of Refs. [5,34], with randomly 
chosen values for the interaction matrix elements that determine the dependence 
of transition energies on the precise population of the single-particle levels in the 
grain, see Subsection 2.1.1. The typical distance between the peaks is of order A/^/g, 
where g is the dimensionless conductance of the metal grain. 

the relation between the splitting of the ground state Kramers degeneracy (in 
the case of a strong spin-orbit coupling and a weak magnetic field) and an 
effective g-tensor. The joint probability distribution of the eigenvalues of the 
tensor was presented. Recent experiments [25] that measured the distribution 
of the (^-tensor in particles of Cu, Ag, and Au, found good agreement with 
many aspects of the theoretical predictions. 

The combined effects of spin-orbit coupling and electron-electron interactions 
were discussed in Subsection 3.2. It was argued that strong spin-orbit cou- 
pling will tend to inhibit the appearance of effects due to electron-electron 
interactions. 

In Subsection 3.3, we reviewed how the peculiar form of the spin-orbit coupling 
for a two-dimensional electron system in a GaAs heterostructure of quantum 
well leads to a strong suppression of spin-orbit effects when the electrons are 
confined in a small quantum dot. We explained how a magnetic field, parallel 



32 



to a quantum dot in a AlGaAs 2DEG, enhances the weak spin-orbit effects in 
these dots. This observation may be used to tune the strength of spin-orbit 
coupling in quantum dots and may explain recent observations on fluctuations 
in the conductance through such dots. 

We also discussed, in Subsection 4 possible explanations, based on non equi- 
librium phenomena, and an almost degenerate ground states due to spin-orbit 
coupling and electron-electron interaction, for the observations [4] of a multi- 
plet splitting of the lowest resonance in the tunneling conductance through a 
gold nanoparticle. 

A few recent developments in the young field of spin and interaction effects in 
small quantum systems have been examined in this article. At present, there 
is no quantitative theory that can explain all the experimental observations 
in this area. Current theories describe well several aspects of small dots (for 
example the (/-tensor eigenvalue distribution), but in others aspects, such as 
Coulomb blockade peak spacing, the agreement between theory and experi- 
ment is far from being satisfactory. 
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A Derivation of the effective Hamiltonian [Eq. (1)] from the toy 
model with contact interaction [Eq. (2)] 

In order to analyze the ground state energy for the toy model Hamiltonian (2), 
the electron-electron interaction is separated into mean and fluctuations, 

uMJ2 n m] n ml , = uM ^ ( (n m] )n ml + n m] (n ml )) 

m m 

+ uM^25n m ^5n ml 

m 
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-«M^(n mt )(n ml ), 

m 



(A.l) 



where 5n ms = n rns — (n ms ) and the average occupation (n m ^) = (n m ^) is 
calculated using the self-consistent (Hartree-Fock) Hamiltonian 

^hf = j2 c t s no(n,m)c miS + uMj2((n m ^n ml + n m ^(n ml )). (A.2) 

n,m,s m 



For technical convenience, we define the occupancies (n ms ) in a reference state 
with 2N — 2K electrons and zero spin, where K is a number of order unity, 
chosen such that all relevant particle-hole excitations have energy less than 
KA. For a disordered metal grain, we may assume that the eigenvectors and 
eigenvalues e^ F of H HF are distributed like those of a random matrix (with 
the possible exception of the spacing of two eigenvalues closest to the Fermi 
level, see Refs. [3,18]; a reference state with 2N — 2K electrons, rather than 
with 2N electrons is chosen, to ensure that only eigenvalues of H RF above the 
Fermi level are needed). The last term on the r.h.s. of Eq. (A.l) is a constant 
shift of the energy and can be omitted. 

We then construct a state with 2N (or 2N + 1) electrons by adding 2K (or 
2K + 1) electrons to the reference state, and find an effective Hamiltonian for 
low-lying particle-hole excitations using the remaining third term on the r.h.s. 
of Eq. (A.l) as a perturbation. The derivation of this effective Hamiltonian 
proceeds along the lines sketched in Sec. 2.2 and Ref. [2]. In the limit M — > oo, 
the effective Hamiltonian has the form (1), where to lowest nontrivial order in 
u one has e^ = e KF and — J s = J c = u. 

When studying the effective interaction amplitudes perturbatively in u, one 
finds that virtual particle-hole excitations involving states with energies e far 
away from E F (A < e — E F < E Th ) contribute to the 0{u 2 ) term and to higher 
order terms. This leads to an effective renormaliziation of the interaction con- 
stants J s and J c , and of the spacing between the levels e^. For example, 



E F 



J s (u) = u + u 2 — 
e»(u) = e™-u 2 ef 



dede 



e + e' 



— I dede 
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e + e' 



+ 0(u 3 ), (A.3) 
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+ 0(u 3 ), 



(A.4) 
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where u(e) is the density of states for the Hamiltonian 7i HF , and J c is renor- 
malized to zero (see appendix C). 

Beyond the first order in the interaction strength u, the symmetry of the effec- 
tive Hamiltonian (1), and of the result (4) that was derived from it, differs from 
that of the equivalent expression in Ref. [3], which was obtained from the toy 
model (2) using the selfconsistent Hartree-Fock approximation. The reason of 
this difference is that, in higher orders of perturbation theory, the selfconsis- 
tent Hartree-Fock approximation neglects certain contributions to the ground 
state energy. (For example, the first correction to J s is of second order in u 
[second term in Eq. (A. 3)], and not, as in the Hartree-Fock approximation, of 
third order [3].) 

When all contributions are taken into account, the symmetry of Eq. (4) and 
the form of the effective Hamiltonian (1) is preserved to all orders in u. 



B The relation between the parameter A = — J s /A and r 



Landau- Fermi-liquid theory expresses various properties of the system in terms 
of the coefficients f pap a'- Using the notations of Ref. [35] we find: 

H = IHL = l + EL. Xp 1 
C b m b d ' XPb m b l + F$ ' 



where C is the specific heat, xp the Pauli susceptibility. The quantities with 
the suffix — b include band effects, but do not include electron-electron in- 
teraction corrections. The latter are encompassed in the F-coefncients of the 
Landau Fermi liquid theory. The letter d denotes the dimension of the system. 
For ballistic systems there are no anomalous renormalizations of the Fermi- 
Liquid coefficients and we have: 

A = - J./A = -F a = 1 - = 1 - (B.2) 

XPb m b xp C b 



Thus, in principle the ratio of the specific heat and the susceptibility gives the 
desired interaction parameter J s . 



B.l Three dimensions 



There are various ways to calculate theoretically the Landau F parameters, 
using different approximations for the electron-electron interaction. They var- 
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ied from a simple static RPA approximation to more complicated approaches 
like density functional theory. Ref. [36] reviews the subject. 



The relevant parameter to describe the strength of the interaction effects is r s , 
the ratio of the typical potential energy to the kinetic energy of electrons. In 
metals it is defined by: 

47r , „ 1 e 2 mb _ . 



where n is the density of electrons and ob is the Bohr radius in the metal. 
Some values for r s in metals are giving in Ref. [37] (page 74) and in Ref. [38]. 

For < r s < 5 the effective mass[37] ranges between 0.96 < m* /m b < 1.06 
where for small r s (< 3), m* < m&, and for larger r s (> 3), m* > m^. (See 
table VII on page 103 in [37].) 

Using the approximation of Rice [39] for the effective mass and for the sus- 
ceptibility one can roughly approximate 

A R ice(r s ) ~ (3 + r s )/25, for 1< r s < 5. (B.4) 



Another approximation for the susceptibility is given in [36] see page 256. 
Assuming that the effective mass is renormalized as in the Rice approach (i.e. 
not very significant renormalization) we find that: 

Asi ng wi(r s ) ~ (2 + r s )/16, for 1< r s < 5. (B.5) 



For small r s ~ 1 the difference between the estimations is only 15 % while for 
r s ~ 5 it close to 30 %. The second estimate reproduces quite well experimental 
measurements of A in a wide range of metals. The parameter A is determined 
by various experimental methods such as electron spin resonance, spin wave, 
Knight shift and total susceptibility. [See p. 256 of Ref. [36] and reference 
therein for further details.] 

Using Eqs. (B.4) and (B.5) we can estimate the parameter J s in different 
materials; however the estimates are rough and should be taken with a grain 
of salt. Typical values for metallic elements are given in the table. 
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Li Na K Rb Cs Cu Ag Au 


\t>- (r ) 
A RiceV s) 

Asingwi^s) 


3.25 3.93 4.86 5.20 5.62 2.67 3.02 3.01 
25 28 31 33 34 23 24 24 
0.33 0.37 0.43 0.45 0.48 0.29 0.31 0.31 




Be Mg Ca Sr Ba Nb Fe Pb 


ARicc(?'s) 
Asingwi( ? *s) 


1.87 2.66 3.27 3.57 3.71 3.07 2.12 2.30 
0.19 0.23 0.25 0.26 0.27 0.24 0.20 0.21 
0.24 0.29 0.33 0.35 0.36 0.32 0.26 0.27 



Table B.l Estimations for A = —J s /A in selected metals. 



B.2 Two dimensions 



Most of the calculations for the Landau- Fermi-liquid parameters in two dimen- 
sional systems were performed for Silicon MOSFET. For a review see Ref. [36] 
(especially page 257) and Ref. [40] (pages 454 and 468). 

We note that as we sweep an external magnetic field, perpendicular to the 
sample area, the spin susceptibility oscillates because the difference in the 
occupations of Landau levels with spins up and down. This effect makes the 
comparison between theory and experiment complicated. We will not be in- 
terested in such anomalously large exchange enhancement. 

In case of silicon MOSFET we should include also the valley degeneracy, and 
the difference in the dielectric constants of Si and SiO that causes the dielectric 
function be space dependent. The screening from the metallic electrodes may 
influence the results as well. 

For GaAs/AlGaAs heterostructures the first two complications are absent. It 
appears that due to the absence of valley degeneracy in GaAs/AlGaAs the 
parameter J s should be larger than in the case of the Si MOSFET. Therefore, 
GA/AlGaAs might be more appropriate to the study of spin configurations 
and there dependance on interaction constants. 
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A static random phase approximation for GaAs gives [41] 



SarccoshCl/x) for < j 
n VI— x 2 — 

*EE^m f or x> i 



(B.6) 



where in two dimensions 

e 2 m b 1 5.45 * 10 5 



h 2 ^wn e47re J n (cm 2 ) 



(B.7) 



In the last equality we take e = 12.9, m b = 0.067m e and m e is the free electron 
mass. 

It can be verified that 

G(x) ™ 1/2 and G(x) ^ (x/tt) log(2/x). 

The factor 1/2 for large x is due to the spin degeneracy and appears because 
both spins participate in the screening in the RPA approximation. (In case of 
a MOSFET the factor 1/2 is substituted by 1/4 due to the valley degeneracy.) 

The same static RPA approximation gives for the effective mass: 

m b /m* = 1 - {V2/n)r s + rf/2 + (1 - r 2 )G(r s /V2). (B.8) 



Numerically, in this approximation 0.95 < m*/mb < 1. In other words within 
the static RPA approximation the mass renormalization is not very significant. 

Using this approximation we plot — J s as a function of r s and n in Fig. B.l. 
For example, X(n = 0.7 * 10 u cm~ 2 ) ^ A(r s ^ 2) ~ 0.34, 



C Renormalization of the interaction in the Cooper channel 



In Sec. 2.1 we described how to integrate out the interaction between electrons 
at high frequency in the RG sense. The interaction in the Cooper channel 
deserves a special consideration, as we will see below it reduces substantially 
when the temperature decreases. 

To see how it works in practice we look at the Dyson equation for the ver- 
tex part of the interaction in the Cooper channel (for a precise definition of 
the vertex part see Ref. [42] Sec. 33.3). Since the divergencies in the Cooper 
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Fig. B.l. A as a function of the ratio of the typical potential energy to the kinetic 
energy of electrons r s , and as a function of the electron density n, for GaAs/AlGaAs 
heterostructures in the static RPA approximation. 

channel are logarithmic we can write the Dyson equation, for T larger than 
the inverse of the elastic mean free time r, in a RG form[16]: 

dJ c (T)/dl = -J 2 c (l)/A, l = log(E F /T), T>l/r. (C.l) 



Integration of this equation, from E-p to T gives 

Jc(T) = i + (Jc(^) C /A)iog(^ F /r)- (c ' 2) 



This logarithmic suppression of the interaction in the Cooper channel was first 
discussed in Ref. [43] and is known as the Tolmachev- Anderson- Morel log or 
pseudo-electron-potential log. 

In quasi-two-dimensional samples, for T < 1/r the RG equation (C.l) is 
modified to [16]: 

dJ c (T)/dl = A/(<7tt) - J C 2 (T)/A, 1/r > T > E Th . (C.3) 



We have neglected here the effects of the diffusive motion on the other chan- 
nels. The presence of the term l/(irg) slows down the logarithmic reduction in 
the Cooper channel, and physically describes the enhancement of the interac- 
tion between the electrons in the Cooper channel due to their diffusive motion. 
In case of quasi-one-dimensional systems the full Dyson equation should be 
solved [44]. 

Finally, the process of integration of the motion at high frequencies reaches the 
Thouless energy, and we find the effective Hamiltonian (1). We will analyze 
the Cooper channel for energies below the Thouless energy using the contact 
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model (2). In principle, the behavior of the Cooper channel can be solved 
exactly by the method of Richardson [45]. But, to understand qualitatively 
the reduction of the interaction in the Cooper channel, at energy below 
it is sufficient to solve the Dyson equation for the interaction matrix element 
(aa\ u \aa) in the Cooper channel. Formally this equation is: 



(aa\u\aa) = \aa\u \aa) — ; ; , (C.4) 



{aa\ 


u° 


\vv) (uu 


u 


aa) 









with u° = uMc\^c\ A c nA c n ^ nm , u = cl^ nl u{n,m)c mi c mV the operator 

Waid) = ^«=i ^ a T(l)( n ) c i t(|) wnere ^e functions M (n) are real eigenfunc- 
tions of a random matrix with real elements, n runs on the sites of the ran- 
dom systems, M is the total number of sites, and \aa) = V'lj.V'if |0)- Using 
the anti-commutation relations of c n >s operators we find an equation for the 
unknown amplitudes u(n,m) 



J2<Pl{n)u{n,m)<f)l{m) = 

nm 

E#(n)(uM* nm )#M + (C.5) 

nm 

^(n)( M M < y nt )^(Q^(fc) M (fc,m)^(m) 
e — e 

n,m,l,k,Vy£a a v 

Now with the relation (<f) 2 u {l)<i) 2 u {k)) = 1/M 2 (1 + 25 ik ) (that is valid in the 
limit M — > oo) we find, comparing the elements in the series that sum over n 
and m 

u(n, m) = uM5 nm — [ m ) + 2w(n, m) J , (C.6) 



where the logarithmic factor appears from the summation over the energies v. 
The term oc 2u (n, m) on the left hand side can be neglected as it is small [by a 
factor (log M)/M] compared to the one on the right hand side. A summation 
over n gives now 

V- (i \ uM 
> u{L,m) = — 

i 1 + AlogM 

substituting in (C.6) and taking the limit M > 1 we find: 

M5 nm + XlogM{M5 nm -l) 

u(n,m) = — , C.7 

V ' ; 1 + AlogM ' v ; 
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where A = u/A and A is the average level spacing in the dot. Hence: 

1 u 
(aa\ u \aa) = — ^(1 + 25 nm )u{n, m) = 2u + — . (C.8) 

m,n 1 + A log 1V1 

The logarithmic factor reduces the term that is oc J2 u(n, m) and does 
not involve the 5 nm factor. This term correspond to the Cooper channel as it 
involves contraction of two wave functions, associated with two creation (or 
two annihilation) operators. Thus, we find that from the Thouless energy up 
to energies of the order of the level spacing, similarly to the pseudo-electron- 
potential log [see after Eq. (C.2)], there is an additional logarithmic suppres- 
sion of the interaction in the Cooper channel. For that reason we took J c = in 
the analysis of the ground state spin distributions in Sec. 2.2 and the Coulomb 
blockade peak spacing in Sec. 2.3. 
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